### R Code for Study 1, MTurk Data ###

# Load Study1_MTurkData.csv #
ExpData <- read.table(file.choose(), header = TRUE, sep = ",")

# Fracking Attitude, Agree/Disagree#
ExpData$FrAttitude <- ifelse(!is.na(ExpData$Q8), ExpData$Q8, ExpData$Q9)
table(ExpData$FrAttitude)

ExpData$FrackAgree <- ifelse(ExpData$FrAttitude < 4, 1, 0)

# Treatment coded 1 if treatment (episodic), 0 if control (thematic)
ExpData$FrTreat <- ifelse(!is.na(ExpData$Q9), 1, 0)
table(ExpData$FrTreat)

# Extremity coded high numbers if extreme #
ExpData$FrExtremity <- abs(4 - ExpData$FrAttitude)
table(ExpData$FrExtremity)
ExpData$FrExtremity <- ExpData$FrExtremity/3
table(ExpData$FrExtremity)

# Importance coded high numbers if important
ExpData$FrImportance <- 5 - ExpData$Q10
table(ExpData$FrImportance)
ExpData$FrImportance <- ExpData$FrImportance/4
table(ExpData$FrImportance)

# Personal Salience coded high numbers if decision affects you 
ExpData$FrPersonalSalience <- 5 - ExpData$Q11
table(ExpData$FrPersonalSalience)
ExpData$FrPersonalSalience <- ExpData$FrPersonalSalience/4
table(ExpData$FrPersonalSalience)

# Moral Conviction 1 -- High numbers indicate moralized
ExpData$FrMoral1 <- 5 - ExpData$Q12

# Moral Conviction 2 -- High numbers indicate moralized
ExpData$FrMoral2 <- 5 - ExpData$Q13

# Full Moral, High numbers means high moral
ExpData$FrMoralFull <- ExpData$FrMoral1 + ExpData$FrMoral2
table(ExpData$FrMoralFull)
ExpData$FrMoralFull <- ExpData$FrMoralFull/8
table(ExpData$FrMoralFull)

# Correlations 
cor(ExpData$FrMoralFull, ExpData$FrPersonalSalience, use = "pairwise.complete.obs")
cor(ExpData$FrMoralFull, ExpData$FrImportance, use = "pairwise.complete.obs")
cor(ExpData$FrMoralFull, ExpData$FrExtremity, use = "pairwise.complete.obs")

# Create Emo Vars 
ExpData$Anger <- ExpData$Q29_1
ExpData$Disgust <- ExpData$Q29_2
ExpData$Sadness <- ExpData$Q29_3
ExpData$Anxiety <- ExpData$Q29_4

emodat <- c("Anger", "Disgust",
            "Sadness", "Anxiety")

emodata <- ExpData[emodat]

emodata <- na.omit(emodata)

# Coding emo vars, cronbach's alpha
library(psych)
alpha(emodata)

ExpData$AllEmo <- ExpData$Anger + ExpData$Disgust + 
  ExpData$Sadness + ExpData$Anxiety
table(ExpData$AllEmo)
ExpData$AllEmo <- ExpData$AllEmo - 4
table(ExpData$AllEmo)
ExpData$AllEmo <- ExpData$AllEmo/16
table(ExpData$AllEmo)

ExpData$Anger <- (ExpData$Anger - 1)/4
table(ExpData$Anger)
ExpData$Disgust <- (ExpData$Disgust - 1)/4
table(ExpData$Disgust)
ExpData$Sadness <- (ExpData$Sadness - 1)/4
table(ExpData$Sadness)
ExpData$Anxiety <- (ExpData$Anxiety - 1)/4
table(ExpData$Anxiety)

# Coding SUR Mods 
library(systemfit)

e1 <- FrMoralFull ~ FrTreat
e2 <- FrPersonalSalience ~ FrTreat
e3 <- FrImportance ~ FrTreat
e4 <- FrExtremity ~ FrTreat

SURVars <- c("FrMoralFull", "FrImportance", "FrPersonalSalience",
             "FrExtremity", "FrTreat")

ExpData1 <- ExpData[SURVars]
ExpData1 <- na.omit(ExpData1)

fitsur <- systemfit(list(Moral = e1, PerSal = e2, Imp = e3, Ext = e4), 
                    data=ExpData1)
summary(fitsur)

library(car)
rest1 <- "Moral_FrTreat - PerSal_FrTreat"
linearHypothesis(fitsur, rest1, test = "F")

rest2 <- "Moral_FrTreat - Imp_FrTreat"
linearHypothesis(fitsur, rest2, test = "F")

rest3 <- "Moral_FrTreat - Ext_FrTreat"
linearHypothesis(fitsur, rest3, test = "F")

### Do Treatment Effects Vary By Agree/Disagree (Supp Material)

e5 <- FrMoralFull ~ FrTreat*FrackAgree
e6 <- FrPersonalSalience ~ FrTreat*FrackAgree
e7 <- FrImportance ~ FrTreat*FrackAgree
e8 <- FrExtremity ~ FrTreat*FrackAgree

SURVars2 <- c("FrMoralFull", "FrImportance", "FrPersonalSalience",
              "FrExtremity", "FrTreat", "FrackAgree")

ExpData2 <- ExpData[SURVars2]
ExpData3 <- na.omit(ExpData2)

fitsur2 <- systemfit(list(Moral = e5, PerSal = e6, Imp = e7, Ext = e8), 
                     data=ExpData3)
summary(fitsur2)

# Coding for Pol Vars and Demos for Mediation Analyses

ExpData$Ind <- 0
ExpData$Ind[ExpData$Q20 == 3] <- 1
ExpData$Ind[ExpData$Q20 == 4] <- 1

ExpData$Dem <- 0
ExpData$Dem[ExpData$Q20 == 1] <- 1

ExpData$Rep <- 0
ExpData$Rep[ExpData$Q20 == 2] <- 1

table(ExpData$Q24)
ExpData$Ideol <- ExpData$Q24

table(ExpData$Q25)
ExpData$Female <- ifelse(ExpData$Q25 == 2, 1, 0)

table(ExpData$Q26)
ExpData$Age <- ExpData$Q26

MedVars <- c("FrTreat", "FrMoralFull", "AllEmo","FrExtremity",
             "FrImportance", "FrPersonalSalience", "Ind", "Dem",
             "Rep", "Ideol", "Female", "Age", "Anger", "Disgust",
             "Sadness", "Anxiety")

MedData <- ExpData[MedVars]
MedData <- na.omit(MedData)

library(mediation)
set.seed(12345)

# Full Emo Scale -- Moral Conviction #
set.seed(12345)
med.fit.emo1 <- lm(AllEmo ~ FrTreat + Dem + Rep + Ideol + Female + Age, 
                   data = MedData)
out.fit.emo1 <- lm(FrMoralFull ~ AllEmo + FrTreat + Dem + Rep + 
                     Ideol + Age + FrExtremity + FrPersonalSalience +
                     FrImportance, data = MedData)
med.out.emo1 <- mediate(med.fit.emo1, out.fit.emo1, treat = "FrTreat",
                        mediator = "AllEmo", robustSE = TRUE, 
                        sims = 2000)
summary(med.out.emo1)

# Full Emo Scale -- Attitude Extremity #
set.seed(12345)
med.fit.emo6 <- lm(AllEmo ~ FrTreat + Dem + Rep + Ideol + Female + Age, 
                   data = MedData)
out.fit.emo6 <- lm(FrExtremity ~ AllEmo + FrTreat + Dem + Rep + 
                     Ideol + Age + FrMoralFull + FrPersonalSalience +
                     FrImportance, data = MedData)
med.out.emo6 <- mediate(med.fit.emo6, out.fit.emo6, treat = "FrTreat",
                        mediator = "AllEmo", robustSE = TRUE, 
                        sims = 2000)
summary(med.out.emo6)

# Full Emo Scale -- Attitude Importance #
set.seed(12345)
med.fit.emo11 <- lm(AllEmo ~ FrTreat + Dem + Rep + Ideol + Female + Age, 
                    data = MedData)
out.fit.emo11 <- lm(FrImportance ~ AllEmo + FrTreat + Dem + Rep + 
                      Ideol + Age + FrMoralFull + FrExtremity +
                      FrPersonalSalience, data = MedData)
med.out.emo11 <- mediate(med.fit.emo11, out.fit.emo11, treat = "FrTreat",
                         mediator = "AllEmo", robustSE = TRUE, 
                         sims = 2000)
summary(med.out.emo11)

# Full Emo Scale -- Personal Salience
set.seed(12345)
med.fit.emo16 <- lm(AllEmo ~ FrTreat + Dem + Rep + Ideol + Female + Age, 
                    data = MedData)
out.fit.emo16 <- lm(FrPersonalSalience ~ AllEmo + FrTreat + Dem + Rep + 
                      Ideol + Age + FrMoralFull + FrImportance +
                      FrExtremity, data = MedData)
med.out.emo16 <- mediate(med.fit.emo16, out.fit.emo16, treat = "FrTreat",
                         mediator = "AllEmo", robustSE = TRUE, 
                         sims = 2000)
summary(med.out.emo16)

### Multiple Mediators (Supplemental Material) ###

### Moral Conviction ###

# Anger #
set.seed(12345)
m.med_anger_moralconv95 <- multimed(outcome = "FrMoralFull", med.main = "Anger",
                                    med.alt = c("Disgust", "Sadness", "Anxiety"),
                                    treat = "FrTreat",
                                    covariates = c("FrExtremity", "FrImportance",
                                                   "FrPersonalSalience",
                                                   "Dem", "Rep", "Ideol", 
                                                   "Female", "Age"),
                                    data = MedData, sims = 2000,
                                    conf.level = .95)
summary(m.med_anger_moralconv95)

set.seed(12345)
m.med_anger_moralconv99 <- multimed(outcome = "FrMoralFull", med.main = "Anger",
                                    med.alt = c("Disgust", "Sadness", "Anxiety"),
                                    treat = "FrTreat",
                                    covariates = c("FrExtremity", "FrImportance",
                                                   "FrPersonalSalience",
                                                   "Dem", "Rep", "Ideol", 
                                                   "Female", "Age"),
                                    data = MedData, sims = 2000,
                                    conf.level = .99)
summary(m.med_anger_moralconv99)

# Disgust #
set.seed(12345)
m.med_disgust_moralconv95 <- multimed(outcome = "FrMoralFull", med.main = "Disgust",
                                      med.alt = c("Anger", "Sadness", "Anxiety"),
                                      treat = "FrTreat",
                                      covariates = c("FrExtremity", "FrImportance",
                                                     "FrPersonalSalience",
                                                     "Dem", "Rep", "Ideol", 
                                                     "Female", "Age"),
                                      data = MedData, sims = 2000,
                                      conf.level = .95)
summary(m.med_disgust_moralconv95)

set.seed(12345)
m.med_disgust_moralconv99 <- multimed(outcome = "FrMoralFull", med.main = "Disgust",
                                      med.alt = c("Anger", "Sadness", "Anxiety"),
                                      treat = "FrTreat",
                                      covariates = c("FrExtremity", "FrImportance",
                                                     "FrPersonalSalience",
                                                     "Dem", "Rep", "Ideol", 
                                                     "Female", "Age"),
                                      data = MedData, sims = 2000,
                                      conf.level = .99)
summary(m.med_disgust_moralconv99)

# Anxiety #
set.seed(12345)
m.med_anxiety_moralconv95 <- multimed(outcome = "FrMoralFull", med.main = "Anxiety",
                                      med.alt = c("Anger", "Sadness", "Disgust"),
                                      treat = "FrTreat",
                                      covariates = c("FrExtremity", "FrImportance",
                                                     "FrPersonalSalience",
                                                     "Dem", "Rep", "Ideol", 
                                                     "Female", "Age"),
                                      data = MedData, sims = 2000,
                                      conf.level = .95)
summary(m.med_anxiety_moralconv95)

set.seed(12345)
m.med_anxiety_moralconv99 <- multimed(outcome = "FrMoralFull", med.main = "Anxiety",
                                      med.alt = c("Anger", "Sadness", "Disgust"),
                                      treat = "FrTreat",
                                      covariates = c("FrExtremity", "FrImportance",
                                                     "FrPersonalSalience",
                                                     "Dem", "Rep", "Ideol", 
                                                     "Female", "Age"),
                                      data = MedData, sims = 2000,
                                      conf.level = .99)
summary(m.med_anxiety_moralconv99)

# Sadness #
set.seed(12345)
m.med_sadness_moralconv95 <- multimed(outcome = "FrMoralFull", med.main = "Sadness",
                                      med.alt = c("Anger", "Anxiety", "Disgust"),
                                      treat = "FrTreat",
                                      covariates = c("FrExtremity", "FrImportance",
                                                     "FrPersonalSalience",
                                                     "Dem", "Rep", "Ideol", 
                                                     "Female", "Age"),
                                      data = MedData, sims = 2000,
                                      conf.level = .95)
summary(m.med_sadness_moralconv95)

set.seed(12345)
m.med_sadness_moralconv99 <- multimed(outcome = "FrMoralFull", med.main = "Sadness",
                                      med.alt = c("Anger", "Anxiety", "Disgust"),
                                      treat = "FrTreat",
                                      covariates = c("FrExtremity", "FrImportance",
                                                     "FrPersonalSalience",
                                                     "Dem", "Rep", "Ideol", 
                                                     "Female", "Age"),
                                      data = MedData, sims = 2000,
                                      conf.level = .99)
summary(m.med_sadness_moralconv99)


### Personal Salience ###

# Anger #
set.seed(12345)
m.med_anger_persal95 <- multimed(outcome = "FrPersonalSalience", med.main = "Anger",
                                 med.alt = c("Disgust", "Sadness", "Anxiety"),
                                 treat = "FrTreat",
                                 covariates = c("FrExtremity", "FrImportance",
                                                "FrMoralFull",
                                                "Dem", "Rep", "Ideol", 
                                                "Female", "Age"),
                                 data = MedData, sims = 2000,
                                 conf.level = .95)
summary(m.med_anger_persal95)

set.seed(12345)
m.med_anger_persal99 <- multimed(outcome = "FrPersonalSalience", med.main = "Anger",
                                 med.alt = c("Disgust", "Sadness", "Anxiety"),
                                 treat = "FrTreat",
                                 covariates = c("FrExtremity", "FrImportance",
                                                "FrMoralFull",
                                                "Dem", "Rep", "Ideol", 
                                                "Female", "Age"),
                                 data = MedData, sims = 2000,
                                 conf.level = .99)
summary(m.med_anger_persal99)

# Disgust #
set.seed(12345)
m.med_disgust_persal95 <- multimed(outcome = "FrPersonalSalience", med.main = "Disgust",
                                   med.alt = c("Anger", "Sadness", "Anxiety"),
                                   treat = "FrTreat",
                                   covariates = c("FrExtremity", "FrImportance",
                                                  "FrMoralFull",
                                                  "Dem", "Rep", "Ideol", 
                                                  "Female", "Age"),
                                   data = MedData, sims = 2000,
                                   conf.level = .95)
summary(m.med_disgust_persal95)

set.seed(12345)
m.med_disgust_persal99 <- multimed(outcome = "FrPersonalSalience", med.main = "Disgust",
                                   med.alt = c("Anger", "Sadness", "Anxiety"),
                                   treat = "FrTreat",
                                   covariates = c("FrExtremity", "FrImportance",
                                                  "FrMoralFull",
                                                  "Dem", "Rep", "Ideol", 
                                                  "Female", "Age"),
                                   data = MedData, sims = 2000,
                                   conf.level = .99)
summary(m.med_disgust_persal99)

# Anxiety #
set.seed(12345)
m.med_anxiety_persal95 <- multimed(outcome = "FrPersonalSalience", med.main = "Anxiety",
                                   med.alt = c("Anger", "Sadness", "Disgust"),
                                   treat = "FrTreat",
                                   covariates = c("FrExtremity", "FrImportance",
                                                  "FrMoralFull",
                                                  "Dem", "Rep", "Ideol", 
                                                  "Female", "Age"),
                                   data = MedData, sims = 2000,
                                   conf.level = .95)
summary(m.med_anxiety_persal95)

set.seed(12345)
m.med_anxiety_persal99 <- multimed(outcome = "FrPersonalSalience", med.main = "Anxiety",
                                   med.alt = c("Anger", "Sadness", "Disgust"),
                                   treat = "FrTreat",
                                   covariates = c("FrExtremity", "FrImportance",
                                                  "FrMoralFull",
                                                  "Dem", "Rep", "Ideol", 
                                                  "Female", "Age"),
                                   data = MedData, sims = 2000,
                                   conf.level = .99)
summary(m.med_anxiety_persal99)

# Sadness #
set.seed(12345)
m.med_sadness_persal95 <- multimed(outcome = "FrPersonalSalience", med.main = "Sadness",
                                   med.alt = c("Anger", "Anxiety", "Disgust"),
                                   treat = "FrTreat",
                                   covariates = c("FrExtremity", "FrImportance",
                                                  "FrMoralFull",
                                                  "Dem", "Rep", "Ideol", 
                                                  "Female", "Age"),
                                   data = MedData, sims = 2000,
                                   conf.level = .95)
summary(m.med_sadness_persal95)

set.seed(12345)
m.med_sadness_persal99 <- multimed(outcome = "FrPersonalSalience", med.main = "Sadness",
                                   med.alt = c("Anger", "Anxiety", "Disgust"),
                                   treat = "FrTreat",
                                   covariates = c("FrExtremity", "FrImportance",
                                                  "FrMoralFull",
                                                  "Dem", "Rep", "Ideol", 
                                                  "Female", "Age"),
                                   data = MedData, sims = 2000,
                                   conf.level = .99)
summary(m.med_sadness_persal99)


### Attitude Importance ###

# Anger #
set.seed(12345)
m.med_anger_imp95 <- multimed(outcome = "FrImportance", med.main = "Anger",
                              med.alt = c("Disgust", "Sadness", "Anxiety"),
                              treat = "FrTreat",
                              covariates = c("FrExtremity", "FrPersonalSalience",
                                             "FrMoralFull",
                                             "Dem", "Rep", "Ideol", 
                                             "Female", "Age"),
                              data = MedData, sims = 2000,
                              conf.level = .95)
summary(m.med_anger_imp95)

set.seed(12345)
m.med_anger_imp99 <- multimed(outcome = "FrImportance", med.main = "Anger",
                              med.alt = c("Disgust", "Sadness", "Anxiety"),
                              treat = "FrTreat",
                              covariates = c("FrExtremity", "FrPersonalSalience",
                                             "FrMoralFull",
                                             "Dem", "Rep", "Ideol", 
                                             "Female", "Age"),
                              data = MedData, sims = 2000,
                              conf.level = .99)
summary(m.med_anger_imp99)

# Disgust #
set.seed(12345)
m.med_disgust_imp95 <- multimed(outcome = "FrImportance", med.main = "Disgust",
                                med.alt = c("Anger", "Sadness", "Anxiety"),
                                treat = "FrTreat",
                                covariates = c("FrExtremity", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .95)
summary(m.med_disgust_imp95)

set.seed(12345)
m.med_disgust_imp99 <- multimed(outcome = "FrImportance", med.main = "Disgust",
                                med.alt = c("Anger", "Sadness", "Anxiety"),
                                treat = "FrTreat",
                                covariates = c("FrExtremity", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .99)
summary(m.med_disgust_imp99)

# Anxiety #
set.seed(12345)
m.med_anxiety_imp95 <- multimed(outcome = "FrImportance", med.main = "Anxiety",
                                med.alt = c("Anger", "Sadness", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrExtremity", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .95)
summary(m.med_anxiety_imp95)

set.seed(12345)
m.med_anxiety_imp99 <- multimed(outcome = "FrImportance", med.main = "Anxiety",
                                med.alt = c("Anger", "Sadness", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrExtremity", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .99)
summary(m.med_anxiety_imp99)

# Sadness #
set.seed(12345)
m.med_sadness_imp95 <- multimed(outcome = "FrImportance", med.main = "Sadness",
                                med.alt = c("Anger", "Anxiety", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrExtremity", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .95)
summary(m.med_sadness_imp95)

set.seed(12345)
m.med_sadness_imp99 <- multimed(outcome = "FrImportance", med.main = "Sadness",
                                med.alt = c("Anger", "Anxiety", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrExtremity", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .99)
summary(m.med_sadness_imp99)


# Attitude Extremity #

# Anger #
set.seed(12345)
m.med_anger_ext95 <- multimed(outcome = "FrExtremity", med.main = "Anger",
                              med.alt = c("Disgust", "Sadness", "Anxiety"),
                              treat = "FrTreat",
                              covariates = c("FrImportance", "FrPersonalSalience",
                                             "FrMoralFull",
                                             "Dem", "Rep", "Ideol", 
                                             "Female", "Age"),
                              data = MedData, sims = 2000,
                              conf.level = .95)
summary(m.med_anger_ext95)

set.seed(12345)
m.med_anger_ext99 <- multimed(outcome = "FrExtremity", med.main = "Anger",
                              med.alt = c("Disgust", "Sadness", "Anxiety"),
                              treat = "FrTreat",
                              covariates = c("FrImportance", "FrPersonalSalience",
                                             "FrMoralFull",
                                             "Dem", "Rep", "Ideol", 
                                             "Female", "Age"),
                              data = MedData, sims = 2000,
                              conf.level = .99)
summary(m.med_anger_ext99)

# Disgust #
set.seed(12345)
m.med_disgust_ext95 <- multimed(outcome = "FrExtremity", med.main = "Disgust",
                                med.alt = c("Anger", "Sadness", "Anxiety"),
                                treat = "FrTreat",
                                covariates = c("FrImportance", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .95)
summary(m.med_disgust_ext95)

set.seed(12345)
m.med_disgust_ext99 <- multimed(outcome = "FrExtremity", med.main = "Disgust",
                                med.alt = c("Anger", "Sadness", "Anxiety"),
                                treat = "FrTreat",
                                covariates = c("FrImportance", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .99)
summary(m.med_disgust_ext99)

# Anxiety #
set.seed(12345)
m.med_anxiety_ext95 <- multimed(outcome = "FrExtremity", med.main = "Anxiety",
                                med.alt = c("Anger", "Sadness", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrImportance", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .95)
summary(m.med_anxiety_ext95)

set.seed(12345)
m.med_anxiety_ext99 <- multimed(outcome = "FrExtremity", med.main = "Anxiety",
                                med.alt = c("Anger", "Sadness", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrImportance", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .99)
summary(m.med_anxiety_ext99)

# Sadness #
set.seed(12345)
m.med_sadness_ext95 <- multimed(outcome = "FrExtremity", med.main = "Sadness",
                                med.alt = c("Anger", "Anxiety", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrImportance", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .95)
summary(m.med_sadness_ext95)

set.seed(12345)
m.med_sadness_ext99 <- multimed(outcome = "FrExtremity", med.main = "Sadness",
                                med.alt = c("Anger", "Anxiety", "Disgust"),
                                treat = "FrTreat",
                                covariates = c("FrImportance", "FrPersonalSalience",
                                               "FrMoralFull",
                                               "Dem", "Rep", "Ideol", 
                                               "Female", "Age"),
                                data = MedData, sims = 2000,
                                conf.level = .99)
summary(m.med_sadness_ext99)